Neutrino induced pion production at MiniBooNE and K2K 
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We investigate charged and neutral current neutrino induced incoherent pion production off nuclei 
at MiniBooNE and K2K energies within the GiBUU model. We assume impulse approximation and 
treat the nucleus as a local Fermi gas of nucleons bound in a mean-field potential. In-medium 
spectral functions are also taken into account. The outcome of the initial neutrino nucleon reaction 
undergoes complex hadronic final state interactions. We present results for neutral current tt" and 
charged current vr^ production and compare to MiniBooNE and K2K data. 
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Introduction. A proper understanding of neutrino 
induced pion production is essential for the interpreta- 
tion of current neutrino oscillations experiments since it 
plays an important role in neutrino flavor identification. 
7r° production events in neutral current (NC) reactions 
are a source of background in Ve appearance searches be- 
cause they might be misidentified as charge current (CC) 
(i/e,e~) interactions. Similarly, CC induced tt^'" mesons 
are a background in i/^ disappearance experiments. As 
all of the present oscillations experiments use nuclear tar- 
gets, it is mandatory to consider final state interactions 
(FSl), i.e., pion rescattering, with and without charge 
exchange, and absorption in the nuclear medium. Next 
generation experiments such as T2K require the CC 7r+ 
(NC 77°) cross section to be known to 5% (10%) for the 
resulting error on the oscillation parameters to be com- 
parable to that from statistical uncertainties ^ . 

Neutrino induced pion production on nucleons up to 
energies of about 1.5 GeV is dominated by the excitation 
and subsequent decay of the A(1232) resonance but, de- 
pending on the channel, non-resonant pion production 
is not negligible; at higher energies, heavier resonances 
become increasingly important [3, y, 0, H, 0, 0, S] ■ For 
nucleons bound in a nucleus, the cross sections are mod- 
ified due to Fermi motion, Pauli blocking, mean-field po- 
tentials, and coUisional broadening of the particles 
Praet et al. [lO] employ realistic bound-state wave func- 
tions and find a good agreement to relativistic Fermi gas 
calculations at neutrino energies around 1 GeV. 

On nuclei, pions can be produced either coherently, 
leaving the nucleus intact, or incoherently. While the 
former one has attracted considerable attention, the lit- 
erature on incoherent processes is limited. A full de- 
scription of the pion production process requires a re- 
alistic treatment of the FSI in the nucleus as outlined 
before. The study of semi-inclusive pion production pre- 
sented in Ref. [ll| includes pion absorption and charge 
exchange but does not properly take into account the im- 
portant features of TriVA dynamics, leading to unrealistic 
pion spectra 0, [l3|. Extending the model of Ref. [H, 
Ahmad et al. [145 include, besides pion absorption, also 



elastic and charge exchange rescattering using empirical 
vacuum ttN cross sections. FSI are neglected in Ref. [l3| 
who calculate pion spectra in a plane wave impulse ap- 
proximation. 

Monte Carlo event generators, which are used in the 
simulation of the neutrino experiments, are commonly 
based on the Rein-Sehgal model for pion production on 
the nucleon [H but they differ significantl y in the treat- 
ment of nuclear effects and FSI (cf. Ref. [l5| and refer- 
ences therein). First MiniBooNE and K2K results indi- 
cate a disagreement between the Monte Carlo predictions 
and the actual measurements [lit . 

In this brief report, we apply the Giessen Boltzmann- 
Uehling-Uhlenbeck (GiBUU) model for neutrino induced 
reactions to NC and CC pion production on Carbon and 
Oxygen nuclei at K2K and MiniBooNE energies. Gen- 
eral results were presented in Refs. IJ, |l7|. The GiBUU 
model is based on well-founded theoretical ingredients 
and has been tested in various and very different reac- 
tions, in particular, against electron and photon scatter- 
ing data [1, [l3| ■ 

GiBUU model. In our model, neutrino nucleus scat- 
tering is treated as a two step process. In the initial state 
step, the neutrinos interact with nucleons embedded in 
the nuclear medium as explained in detail in Ref. 0]. 
In the final state step, the outgoing particles of the ini- 
tial reaction are propagated through the nucleus, using a 
hadronic transport approach As we shall see, these 
FSI modify the pion yields considerably. 

We treat the nucleus as a local Fermi gas of nucleons 
bound in a mean-field potential and obtain for the total 
pion production cross section on nuclei 
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k is the four-vector of the the neutrino, p the one of the 
bound nucleon, andp^'^(r) = (37r^p"'P(r))^/'^ denotes the 
local Fermi momentum depending on the nuclear density. 
Mtt is the multiplicity of the final state which is deter- 
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mined by the GiBUU transport simulation described be- 
low, dcrj^t'^ stands for the total cross section on nucleons 
including nuclear medium corrections. 

For neutrino beam energies ranging from 0.5 to 2 GeV, 
pions are produced both through resonance and back- 
ground contributions Even initial quasielastic scat- 
tering (QE) events contribute to the pion production 
cross section through secondary collisions in the FSI. The 
pion production cross section is dominated by the exci- 
tation of the A resonance P33(1232). Additionally, we 
include 12 TV* and A resonances with invariant masses 
less than 2 GeV and also non-resonant pion produc- 
tion (treated as background in our model) which is non- 
negligible in some cases. The vector parts of the reso- 
nance contributions are obtained from a recent MAID 
analysis of electroproduction cross sections, the same 
holds for the background vector contributions. The ax- 
ial couplings for the resonances are obtained from PCAC 
(partial conservation of the axial current) as described 
in Ref. We further use neutrino nucleon data to fit 
the axial form factor of the A resonance as well as the 
non-vector background contribution. We refer the reader 
to Ref. Q where our model for lepton-nucleon scattering 
is described in detail. 

The neutrino nucleon cross sections are modified in the 
nuclear medium. Bound nucleons are treated within a lo- 
cal Thomas-Fermi approximation. They are exposed to 
a mean-field potential depending on density and momen- 
tum. We account for this by evaluating the above cross 
sections with full in- medium kinematics, i.e., hadronic 
tensor, flux and phase-space factors are evaluated with 
in-medium four-vectors. We also take Pauli blocking into 
account. 

Once produced inside the nucleus, the particles prop- 
agate out to the detector. During their propaga- 
tion they undergo FSI which are simulated with the 
coupled-channel semi-classical GiBUU transport model 
[igt.^ Originally developed to describe heavy-ion col- 
lisions, it has been extended to describe the interac- 
tions of pion s, photons, electrons, and neutrinos with 
nuclei 1,111, [2^. 

In the following, we give a brief review on the basic in- 
gredients of our model; for full details we refer the reader 
to Refs. and references therein. It is based on the 
BUU equation which describes the space-time evolution 
of a many-particle system in a mean-field potential. For 
particles of species i, it is given by 

{dt + VpH ■ V, - VrH ■ Vp) Mv,p,t) = 

Icoll[fi, In, f-rr, /a, ■■■], (2) 



where the phase space density fi{r,p, t) depends on time 
t, coordinates r and the four-momentum p. H is the rel- 
ativistic Hamiltonian of a particle of mass M in a scalar 

potential U given hy H = \J[M + t/(r,p)]^ +p2. The 
collision term Icoii accounts for changes (gain and loss) in 
the phase space density due to elastic and inelastic col- 
lisions between the particles, and also to particle decays 
into other hadrons. The BUU equations for all parti- 
cle species are thus coupled through the collision term 
and also through the potentials in H. A coupled-channel 
treatment is required to take into account side-feeding 
into different channels. Baryon-meson two-body interac- 
tions (e.g., ttN ttN) are dominated by resonance con- 
tributions and a small non-resonant background term; 
baryon-baryon cross sections (e.g., NN NN, RN 
NN, RN R'N, NN ttNN) are either fitted to 
data or calculated e.g., in pion exchange models. The 
three-body channels ttNN NN and ANN NNN 
are also included. 

All particles (also resonances) are propagated in mean- 
field potentials according to their BUU equation. Those 
states acquire medium-modified spectral functions (nu- 
cleons and resonances) and are propagated off-shell. 
The medium-modification of the spectral function is 
based both on coUisional broadening and on the mean- 
field potentials. The collisional broadening is cal- 
culated using the low-density approximation Fcoii = 
a{E,p,p')vreip{r), where a{E, p,p') denotes the to- 
tal cross section for the scattering of the outgoing nu- 
cleon/resonance of energy E and momentum p with a 
nucleon of momentum p' in the vacuum with relative ve- 
locity Vi-ci- The collisional broadening is obtained in a 
consistent way from the GiBUU cross sections, and de- 
pends on the particle kinematics as well as on the nuclear 
density. In Ref. 0, [2l[, we have shown that both the 
momentum dependence of the mean-field potential and 
the collisional broadening are necessary to obtain good 
agreement with inclusive electron scattering experimen- 
tal data. We ensure that after leaving the nucleus, vac- 
uum spectral functions are recovered. Finally, the pion 
multiplicity M^^ is determined by counting all asymptotic 
pions in each kinematical bin. 

Summarizing, FSI lead to absorption, charge exchange 
and redistribution of energy and momentum, as well as 
to the production of new particles. Their impact on ncu- 



^ The numerical implementation of the GiBUU model is available 
for download from our website [T9II . 



trino induced pion production is dramatic |13l . Il7l | and, 
therefore, a qualitatively and quantitatively correct treat- 
ment is of great importance. 

Results. In Fig. [T] (a), we show our results for NC 
single- 7r° production on ^^C as a function of the pion 
kinetic energy. We have averaged over the MiniBooNE 
energyflux which peaks at about 0.7 GeV neutrino en- 
ergy [22]. As can be deduced from simple isospin argu- 
ments, in NC reactions the total pion yield is dominated 
by 7r° production, while 7r+ dominate in CC processes. 
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Comparing the dashed with the sohd hne (results with- 
out FSI and spectral function vs. full calculation), one 
finds a considerably change of the spectra. The shape 
is caused by the energy dependence of the pion absorp- 
tion and rescattering cross sections. Pions are mainly 
absorbed via the A resonance, i.e, through ttN A fol- 
lowed by AA^ NN. This explains the reduction in the 
region around T^. = 0.1 — 0.3 GeV. Pion elastic scatter- 
ing ttN ttN reshuffles the pions to lower momenta and 
leads also to charge exchange scattering into the charged 
pion channels. The vast majority of the pions come from 
initial A excitation (dash-dotted line). Their production 
in the rescattering of nucleons is not significant here but 
becomes more important at higher energies and for heav- 
ier nuclei [3]. 

The MiniBooNE experiment has recently measured 
NC single- tt" momentum spectra [2^, however, their data 
are available only as count rates. Notice that the data 
include a small contribution from coherent pion produc- 
tion, i.e., vA v'k'^A, which cannot be described by our 
transport model. A direct and meaningful comparison 
to the MiniBooNE measurement will be possible when 
acceptance corrected cross sections are provided. 

Fig.[T](b) shows the results for NC single- 7r° production 
off ^^O averaged over the K2K ener gy f lux which peaks 
at about 1.2 GeV neutrino energy [2J]. Compared to 
(a), the spectrum is broader and extends to larger 
due to the higher neutrino energy. The reduction in the 
region around — 0.1 — 0.3 GeV is mainly caused by the 
pion absorption via the A resonance (compare dashed 
and solid lines). Again, pion production through initial 
QE scattering is not sizable. 

NC single- 7r° spectra have been measured by the K2K 
collaboration [2J|. Based on their Monte Carlo gener- 
ator, the data are not only corrected for efficiency, but 
some background has already been subtracted, i.e., the 
data include a model dependence. Since K2K, as well 
as MiniBooNE, has not yet provided cross sections but 
only count rates, we cannot yet compare to these mea- 
surements directly. 

In Fig. [21 we give our results for the single- 7r+/QE 
ratio for CC interactions on ^^C. We present different 
scenarios: 

(1) cri^+/aQ^+ after FSI: This denotes our full calcu- 
lation after FSI for the single-7r+ cross section di- 
vided by the "CCQE-like" cross section, i.e., the 
cross section after FSI without any pion leaving 
the nucleus. This definition of "CCQE-like" as all 
the events where no pion is detected, is used, e.g., 
by MiniBooNE. 

(2) <Ti^+ /(TQ7r+ ip after FSI: similar to (1), but we re- 
quire in addition a single-p in the final state for the 
"CCQE-like" cross section. Here, "CCQE-like" de- 
notes all events where a single proton track is visi- 
ble and, at the same time, no pions are detected as 
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FIG. 1: (Color online) (a) NC induced single-Tr" production 
on ^^C as a function of the pion kinetic energy averaged over 
the MiniBooNE flux, (b) same on ^''O averaged over the 
K2K flux. The dashed line shows our calculation without FSI 
or spectral functions, both included in the full calculation 
denoted by the solid line. The dash-dotted line indicates the 
contribution from the A resonance to the full calculation. 



apphed, e.g., by K2K. 

(3) (71^+ after FSI/aqE: Here, the full single-7r+ cross 
section after FSI is divided by the total QE cross 
section. 

(4) ai^+ before FSI/ctqe: Same as (3), but before FSI; 
nuclear corrections like mean fields and Fermi mo- 
tion are included. 

(5) (71^+ /(Tqe in the vacuum: For comparison, we show 
the vacuum cross section on an isoscalar target. 

In (3)-(5) (Tqe defines the situation in which the pri- 
mary neutrino-nucleon interaction is purely quasielastic 
{v^n —>■ fJ-~p) regardless of the fate of the struck nucleon. 

We want to emphasize two issues: First, nuclear cor- 
rections cancel out in the ratio, as long as FSI are not 
considered ((4) vs. (5)). In general, the complexity of 
FSI prevent such cancellations as one can infer from (1) 
and (3) . Only on very specific occasions, FSI effects may 
cancel: (2) with the particular definition of the "CCQE- 
like" cross section (= cro7r+ ip) lies nearly on top of (4) 
and (5). 

We further compare to preliminary MiniBooNE data 
(Fig.HKa)) and to K2K data M (Fig.[l](b)). These 
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FIG. 2: (Color online) Single- tt^/QE cross section ratio for 
CC interactions on ^'^C vs. neutrino energy. Our results are 
the same in both panels. The meaning of the five different 
curves is described in the text, (a) shows preliminary data 
from MiniBooNE [H; (b) recent K2K data t26|]. The curves 
(2), (4), and (5) lie nearly on top of each other. 



data are corrected for FSI using specific Monte Carlo gen- 
erators, i.e., they give the cross sections for bound nucle- 
ons "before FSI" . As this procedure introduces a model 
dependence in the data, a fully consistent comparison 
is not possible. Furthermore, also the reconstruction of 
the neutrino energy out of the observed muon and the 
hadrons is model dependent. Ignoring the model depen- 
dencies, our calculation denoted by (4) should be the one 
to compare with. In the MiniBooNE case, the agreement 
is perfect for energies up to 1 GeV, and still within their 
error bars above 1 GeV (cf. Fig. [5] (a)). Wc also reach 
a good agreement with the K2K data (cf. Fig. [2] (b)). 
The slight underestimate of the pion/quasielastic ratio at 
higher energies could be due to either an underestimate 
of the pion production cross section, an overestimate of 
the QE cross section, or to uncertainties in the "data" ex- 
traction which involves a model dependence (see above) . 

To summarize, in this brief report, we have presented 
results for NC and CC single-7r production for Mini- 
BooNE and K2K energies. In particular, we have in- 
vestigated the effect of final state interactions. Wherever 



possible, we have compared our calculation to recent data 
and we have found good agreement. However, all these 
"data" were readjusted using specific Monte Carlo event 
generator with specific assumptions on the initial neu- 
trino nucleon cross section and the nuclear model. Model 
independent acceptance corrected data are required to 
perform meaningful comparisons with theoretical calcu- 
lations. 
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